Genomic insights into anthropozoonotic tuberculosis in captive sun bears (Helarctos malayanus) and an Asiatic black bear (Ursus thibetanus) in Cambodia.

Contact between humans and wildlife presents a risk for both zoonotic and anthropozoonotic disease transmission. In this study we report the detection of human strains of Mycobacterium tuberculosis in sun bears and an Asiatic black bear in a wildlife rescue centre in Cambodia, confirming for the first time the susceptibility of these bear species to tuberculosis when in close contact with humans. After genotyping revealed two different strains of M. tuberculosis from cases occurring between 2009 and 2019, 100 isolates from 30 sun bear cases, a single Asiatic black bear case, and a human case were subjected to whole genome sequencing. We combined single nucleotide polymorphism analysis and exploration of mixed base calls with epidemiological data to indicate the evolution of each outbreak. Our results confirmed two concurrent yet separate tuberculosis outbreaks and established a likely transmission route in one outbreak where the human case acted as an intermediatory between bear cases. In both outbreaks, we observed high rates of transmission and progression to active disease, suggesting that sun bears are highly susceptible to tuberculosis if exposed under these conditions. Overall, our findings highlight the risk of bi-directional transmission of tuberculosis between humans and captive bears in high human tuberculosis burden regions, with implied considerations for veterinary and public health. We also demonstrate the use of standard genomic approaches to better understand disease outbreaks in captive wildlife settings and to inform control and prevention measures.

Targeted sequencing from cerebrospinal fluid for rapid identification of drug-resistant tuberculous meningitis.

Mortality from tuberculous meningitis (TBM) remains around 30%, with most deaths occurring within 2 months of starting treatment. Mortality from drug-resistant strains is higher still, making early detection of drug resistance (DR) essential. Targeted next-generation sequencing (tNGS) produces high read depths, allowing the detection of DR-associated alleles with low frequencies. We applied Deeplex Myc-TB-a tNGS assay-to cerebrospinal fluid (CSF) samples from 72 adults with microbiologically confirmed TBM and compared its genomic drug susceptibility predictions to a composite reference standard of phenotypic susceptibility testing (pDST) and whole genome sequencing, as well as to clinical outcomes. Deeplex detected Mycobacterium tuberculosis complex DNA in 24/72 (33.3%) CSF samples and generated full DR reports for 22/24 (91.7%). The read depth generated by Deeplex correlated with semi-quantitative results from MTB/RIF Xpert. Alleles with <20% frequency were seen at canonical loci associated with first-line DR. Disregarding these low-frequency alleles, Deeplex had 100% concordance with the composite reference standard for all drugs except pyrazinamide and streptomycin. Three patients had positive CSF cultures after 30 days of treatment; reference tests and Deeplex identified isoniazid resistance in two, and Deeplex alone identified low-frequency rifampin resistance alleles in one. Five patients died, of whom one had pDST-identified pyrazinamide resistance. tNGS on CSF can rapidly and accurately detect drug-resistant TBM, but its application is limited to those with higher bacterial loads. In those with lower bacterial burdens, alternative approaches need to be developed for both diagnosis and resistance detection.

Rifampicin resistant Mycobacterium tuberculosis in Vietnam, 2020-2022.

Objective: We conducted a descriptive analysis of multi-drug resistant tuberculosis (MDR-TB) in Vietnam's two largest cities, Hanoi and Ho Chi Minh city. Methods: All patients with rifampicin resistant tuberculosis were recruited from Hanoi and surrounding provinces between 2020 and 2022. Additional patients were recruited from Ho Chi Minh city over the same time period. Demographic data were recorded from all patients, and samples collected, cultured, whole genome sequenced and analysed for drug resistance mutations. Genomic susceptibility predictions were made on the basis of the World Health Organization's catalogue of mutations in Mycobacterium tuberculosis associated with drug resistance, version 2. Comparisons were made against phenotypic drug susceptibility test results where these were available. Multivariable logistic regression was used to assess risk factors for previous episodes of tuberculosis. Results: 233/265 sequenced isolates were of sufficient quality for analysis, 146 (63 %) from Ho Chi Minh City and 87 (37 %) from Hanoi. 198 (85 %) were lineage 2, 20 (9 %) were lineage 4, and 15 (6 %) were lineage 1. 17/211 (8 %) for whom HIV status was known were infected, and 109/214 (51 %) patients had had a previous episode of tuberculosis. The main risk factor for a previous episode was HIV infection (odds ratio 5.1 (95 % confidence interval 1.3-20.0); p = 0.021). Sensitivity for predicting first-line drug resistance from whole genome sequencing data was over 90 %, with the exception of pyrazinamide (85 %). For moxifloxacin and amikacin it was 50 % or less. Among rifampicin-resistant isolates, prevalence of resistance to each non-first-line drug was < 20 %. Conclusions: Drug resistance among most MDR-TB strains in Vietnam's two largest cities is confined largely to first-line drugs. Living with HIV is the main risk factor among patients with MDR-TB for having had a previous episode of tuberculosis.

Prediction of pyrazinamide resistance in Mycobacterium tuberculosis using structure-based machine-learning approaches.

Background: Pyrazinamide is one of four first-line antibiotics used to treat tuberculosis; however, antibiotic susceptibility testing for pyrazinamide is challenging. Resistance to pyrazinamide is primarily driven by genetic variation in pncA, encoding an enzyme that converts pyrazinamide into its active form. Methods: We curated a dataset of 664 non-redundant, missense amino acid mutations in PncA with associated high-confidence phenotypes from published studies and then trained three different machine-learning models to predict pyrazinamide resistance. All models had access to a range of protein structural-, chemical- and sequence-based features. Results: The best model, a gradient-boosted decision tree, achieved a sensitivity of 80.2% and a specificity of 76.9% on the hold-out test dataset. The clinical performance of the models was then estimated by predicting the binary pyrazinamide resistance phenotype of 4027 samples harbouring 367 unique missense mutations in pncA derived from 24 231 clinical isolates. Conclusions: This work demonstrates how machine learning can enhance the sensitivity/specificity of pyrazinamide resistance prediction in genetics-based clinical microbiology workflows, highlights novel mutations for future biochemical investigation, and is a proof of concept for using this approach in other drugs.

Spatial Analysis of Drug-Susceptible and Multidrug-Resistant Cases of Tuberculosis, Ho Chi Minh City, Vietnam, 2020-2023.

We characterized the spatial distribution of drug-susceptible (DS) and multidrug-resistant (MDR) tuberculosis (TB) cases in Ho Chi Minh City, Vietnam, a major metropolis in southeastern Asia, and explored demographic and socioeconomic factors associated with local TB burden. Hot spots of DS and MDR TB incidence were observed in the central parts of Ho Chi Minh City, and substantial heterogeneity was observed across wards. Positive spatial autocorrelation was observed for both DS TB and MDR TB. Ward-level TB incidence was associated with HIV prevalence and the male proportion of the population. No ward-level demographic and socioeconomic indicators were associated with MDR TB case count relative to total TB case count. Our findings might inform spatially targeted TB control strategies and provide insights for generating hypotheses about the nature of the relationship between DS and MDR TB in Ho Chi Minh City and the wider southeastern region of Asia.

Tuberculosis in Ukrainian War Refugees and Migrants in the Czech Republic and Slovakia: A Molecular Epidemiological Study.

BACKGROUND: The war in Ukraine has led to significant migration to neighboring countries, raising public health concerns. Notable tuberculosis (TB) incidence rates in Ukraine emphasize the immediate requirement to prioritize approaches that interrupt the spread and prevent new infections. METHODS: We conducted a prospective genomic surveillance study to assess migration's impact on TB epidemiology in the Czech Republic and Slovakia. Mycobacterium tuberculosis isolates from Ukrainian war refugees and migrants, collected from September 2021 to December 2022 were analyzed alongside 1574 isolates obtained from Ukraine, the Czech Republic, and Slovakia. RESULTS: Our study revealed alarming results, with historically the highest number of Ukrainian tuberculosis patients detected in the host countries. The increasing number of cases of multidrug-resistant TB, significantly linked with Beijing lineage 2.2.1 (p < 0.0001), also presents substantial obstacles to control endeavors. The genomic analysis identified the three highly related genomic clusters, indicating the recent TB transmission among migrant populations. The largest clusters comprised war refugees diagnosed in the Czech Republic, TB patients from various regions of Ukraine, and incarcerated individuals diagnosed with pulmonary TB specialized facility in the Kharkiv region, Ukraine, pointing to a national transmission sequence that has persisted for over 14 years. CONCLUSIONS: The data showed that most infections were likely the result of reactivation of latent disease or exposure to TB before migration rather than recent transmission occurring within the host country. However, close monitoring, appropriate treatment, careful surveillance, and social support are crucial in mitigating future risks, though there is currently no evidence of local transmission in EU countries.

The effect of M. tuberculosis lineage on clinical phenotype.

Six lineages of Mycobacterium tuberculosis sensu stricto (which excludes M. africanum) are described. Single-country or small observational data suggest differences in clinical phenotype between lineages. We present strain lineage and clinical phenotype data from 12,246 patients from 3 low-incidence and 5 high-incidence countries. We used multivariable logistic regression to explore the effect of lineage on site of disease and on cavities on chest radiography, given pulmonary TB; multivariable multinomial logistic regression to investigate types of extra-pulmonary TB, given lineage; and accelerated failure time and Cox proportional-hazards models to explore the effect of lineage on time to smear and culture-conversion. Mediation analyses quantified the direct effects of lineage on outcomes. Pulmonary disease was more likely among patients with lineage(L) 2, L3 or L4, than L1 (adjusted odds ratio (aOR) 1.79, (95% confidence interval 1.49-2.15), p<0.001; aOR = 1.40(1.09-1.79), p = 0.007; aOR = 2.04(1.65-2.53), p<0.001, respectively). Among patients with pulmonary TB, those with L1 had greater risk of cavities on chest radiography versus those with L2 (aOR = 0.69(0.57-0.83), p<0.001) and L4 strains (aOR = 0.73(0.59-0.90), p = 0.002). L1 strains were more likely to cause osteomyelitis among patients with extra-pulmonary TB, versus L2-4 (p = 0.033, p = 0.008 and p = 0.049 respectively). Patients with L1 strains showed shorter time-to-sputum smear conversion than for L2. Causal mediation analysis showed the effect of lineage in each case was largely direct. The pattern of clinical phenotypes seen with L1 strains differed from modern lineages (L2-4). This has implications for clinical management and could influence clinical trial selection strategies.

FLASH-TB: an Application of Next-Generation CRISPR to Detect Drug Resistant Tuberculosis from Direct Sputum.

Offering patients with tuberculosis (TB) an optimal and timely treatment regimen depends on the rapid detection of Mycobacterium tuberculosis (Mtb) drug resistance from clinical samples. Finding Low Abundance Sequences by Hybridization (FLASH) is a technique that harnesses the efficiency, specificity, and flexibility of the Cas9 enzyme to enrich targeted sequences. Here, we used FLASH to amplify 52 candidate genes probably associated with resistance to first- and second-line drugs in the Mtb reference strain (H37Rv), then detect drug resistance mutations in cultured Mtb isolates, and in sputum samples. 92% of H37Rv reads mapped to Mtb targets, with 97.8% of target regions covered at a depth ≥ 10X. Among cultured isolates, FLASH-TB detected the same 17 drug resistance mutations as whole genome sequencing (WGS) did, but with much greater depth. Among the 16 sputum samples, FLASH-TB increased recovery of Mtb DNA compared with WGS (from 1.4% [IQR 0.5-7.5] to 33% [IQR 4.6-66.3]) and average depth reads of targets (from 6.3 [IQR 3.8-10.5] to 1991 [IQR 254.4-3623.7]). FLASH-TB identified Mtb complex in all 16 samples based on IS1081 and IS6110 copies. Drug resistance predictions for 15/16 (93.7%) clinical samples were highly concordant with phenotypic DST for isoniazid, rifampicin, amikacin, and kanamycin [15/15 (100%)], ethambutol [12/15 (80%)] and moxifloxacin [14/15 (93.3%)]. These results highlighted the potential of FLASH-TB for detecting Mtb drug resistance from sputum samples.

The effect of M. tuberculosis lineage on clinical phenotype.

Eight lineages of Mycobacterium tuberculosis sensu stricto are described. Single-country or small observational data suggest differences in clinical phenotype between lineages. We present strain lineage and clinical phenotype data from 12,246 patients from 3 low-incidence and 5 high-incidence countries. We used multivariable logistic regression to explore the effect of lineage on site of disease and on cavities on chest radiography, given pulmonary TB; multivariable multinomial logistic regression to investigate types of extra-pulmonary TB, given lineage; and accelerated failure time and Cox proportional-hazards models to explore the effect of lineage on time to smear and culture-conversion. Mediation analyses quantified the direct effects of lineage on outcomes. Pulmonary disease was more likely among patients with lineage(L) 2, L3 or L4, than L1 (adjusted odds ratio (aOR) 1.79, (95% confidence interval 1.49-2.15), p<0.001; aOR=1.40(1.09-1.79), p=0.007; aOR=2.04(1.65-2.53), p<0.001, respectively). Among patients with pulmonary TB, those with L1 had greater risk of cavities on chest radiography versus those with L2 (aOR=0.69(0.57-0.83), p<0.001) and L4 strains (aOR=0.73(0.59-0.90), p=0.002). L1 strains were more likely to cause osteomyelitis among patients with extra-pulmonary TB, versus L2-4 (p=0.033, p=0.008 and p=0.049 respectively). Patients with L1 strains showed shorter time-to-sputum smear conversion than for L2. Causal mediation analysis showed the effect of lineage in each case was largely direct. The pattern of clinical phenotypes seen with L1 strains differed from modern lineages (L2-4). This has implications for clinical management and could influence clinical trial selection strategies.

Updating the WHO target product profile for next-generation Mycobacterium tuberculosis drug susceptibility testing at peripheral centres.

There were approximately 10 million tuberculosis (TB) cases in 2020, of which 500,000 were drug-resistant. Only one third of drug-resistant TB cases were diagnosed and enrolled on appropriate treatment, an issue partly driven by a lack of rapid, accurate drug-susceptibility testing (DST) tools deployable in peripheral settings. In 2014, World Health Organization (WHO) published target product profiles (TPPs) which detailed minimal and optimal criteria to address high-priority TB diagnostic needs, including DST. Since then, the TB community's needs have evolved; new treatment regimens, changes in TB definitions, further emergence of drug resistance, technological advances, and changing end-users requirements have necessitated an update. The DST TPP's revision was therefore undertaken by WHO with the Stop TB Partnership New Diagnostics Working Group. We describe the process of updating the TPP for next-generation TB DST for use at peripheral centres, highlight key updates, and discuss guidance regarding technical and operational specifications.

Genomic Sequencing from Sputum for Tuberculosis Disease Diagnosis, Lineage Determination, and Drug Susceptibility Prediction.

Universal access to drug susceptibility testing for newly diagnosed tuberculosis patients is recommended. Access to culture-based diagnostics remains limited, and targeted molecular assays are vulnerable to emerging resistance mutations. Improved protocols for direct-from-sputum Mycobacterium tuberculosis sequencing would accelerate access to comprehensive drug susceptibility testing and molecular typing. We assessed a thermo-protection buffer-based direct-from-sample M. tuberculosis whole-genome sequencing protocol. We prospectively analyzed 60 acid-fast bacilli smear-positive clinical sputum samples in India and Madagascar. A diversity of semiquantitative smear positivity-level samples were included. Sequencing was performed using Illumina and MinION (monoplex and multiplex) technologies. We measured the impact of bacterial inoculum and sequencing platforms on genomic read depth, drug susceptibility prediction performance, and typing accuracy. M. tuberculosis was identified by direct sputum sequencing in 45/51 samples using Illumina, 34/38 were identified using MinION-monoplex sequencing, and 20/24 were identified using MinION-multiplex sequencing. The fraction of M. tuberculosis reads from MinION sequencing was lower than from Illumina, but monoplexing grade 3+ samples on MinION produced higher read depth than Illumina (P < 0.05) and MinION multiplexing (P < 0.01). No significant differences in sensitivity and specificity of drug susceptibility predictions were seen across sequencing modalities or within each technology when stratified by smear grade. Illumina sequencing from sputum accurately identified 1/8 (rifampin) and 6/12 (isoniazid) resistant samples, compared to 2/3 (rifampin) and 3/6 (isoniazid) accurately identified with Nanopore monoplex. Lineage agreement levels between direct and culture-based sequencing were 85% (MinION-monoplex), 88% (Illumina), and 100% (MinION-multiplex). M. tuberculosis direct-from-sample whole-genome sequencing remains challenging. Improved and affordable sample treatment protocols are needed prior to clinical deployment.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia.

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

Evaluation of Nanopore sequencing for Mycobacterium tuberculosis drug susceptibility testing and outbreak investigation: a genomic analysis.

BACKGROUND: Mycobacterium tuberculosis whole-genome sequencing (WGS) has been widely used for genotypic drug susceptibility testing (DST) and outbreak investigation. For both applications, Illumina technology is used by most public health laboratories; however, Nanopore technology developed by Oxford Nanopore Technologies has not been thoroughly evaluated. The aim of this study was to determine whether Nanopore sequencing data can provide equivalent information to Illumina for transmission clustering and genotypic DST for M tuberculosis. METHODS: In this genomic analysis, we analysed 151 M tuberculosis isolates from Madagascar, South Africa, and England, which were collected between 2011 and 2018, using phenotypic DST and matched Illumina and Nanopore data. Illumina sequencing was done with the MiSeq, HiSeq 2500, or NextSeq500 platforms and Nanopore sequencing was done on the MinION or GridION platforms. Using highly reliable PacBio sequencing assemblies and pairwise distance correlation between Nanopore and Illumina data, we optimise Nanopore variant filters for detecting single-nucleotide polymorphisms (SNPs; using BCFtools software). We then used those SNPs to compare transmission clusters identified by Nanopore with the currently used UK Health Security Agency Illumina pipeline (COMPASS). We compared Illumina and Nanopore WGS-based DST predictions using the Mykrobe software and mutation catalogue. FINDINGS: The Nanopore BCFtools pipeline identified SNPs with a median precision of 99·3% (IQR 99·1-99·6) and recall of 90·2% (88·1-94·2) compared with a precision of 99·6% (99·4-99·7) and recall of 91·9% (87·6-98·6) using the Illumina COMPASS pipeline. Using a threshold of 12 SNPs for putative transmission clusters, Illumina identified 98 isolates as unrelated and 53 as belonging to 19 distinct clusters (size range 2-7). Nanopore reproduced 15 out of 19 clusters perfectly; two clusters were merged into one cluster, one cluster had a single sample missing, and one cluster had an additional sample adjoined. Illumina-based clusters were also closely replicated using a five SNP threshold and clustering accuracy was maintained using mixed Illumina and Nanopore datasets. Genotyping resistance variants with Nanopore was highly concordant with Illumina, having zero discordant SNPs across more than 3000 SNPs and four insertions or deletions (indels), across 60 000 indels. INTERPRETATION: Illumina and Nanopore technologies can be used independently or together by public health laboratories performing M tuberculosis genotypic DST and outbreak investigations. As a result, clinical and public health institutions making decisions on which sequencing technology to adopt for tuberculosis can base the choice on cost (which varies by country), batching, and turnaround time. FUNDING: Academy for Medical Sciences, Oxford Wellcome Institutional Strategic Support Fund, and the Swiss South Africa Joint Research Award (Swiss National Science Foundation and South African National Research Foundation).

In silico evaluation of WHO-endorsed molecular methods to detect drug resistant tuberculosis.

Universal drug susceptibility testing (DST) for tuberculosis is a major goal of the END TB strategy. PCR-based molecular diagnostic tests have been instrumental in increasing DST globally and several assays have now been endorsed by the World Health Organization (WHO) for use in the diagnosis of drug resistance. These endorsed assays, however, each interrogate a limited number of mutations associated with resistance, potentially limiting their sensitivity compared to sequencing-based methods. We applied an in silico method to compare the sensitivity and specificity of WHO-endorsed molecular based diagnostics to the mutation set identified by the WHO mutations catalogue using phenotypic DST as the reference. We found that, in silico, the mutation sets used by probe-based molecular diagnostic tests to identify rifampicin, isoniazid, pyrazinamide, levofloxacin, moxifloxacin, amikacin, capreomycin and kanamycin resistance produced similar sensitivities and specificities to the WHO mutation catalogue. PCR-based diagnostic tests were most sensitive for drugs where mechanisms of resistance are well established and localised to small genetic regions or a few prevalent mutations. Approaches using sequencing technologies can provide advantages for drugs where our knowledge of resistance is limited, or where complex resistance signatures exist.

A crowd of BashTheBug volunteers reproducibly and accurately measure the minimum inhibitory concentrations of 13 antitubercular drugs from photographs of 96-well broth microdilution plates.

Tuberculosis is a respiratory disease that is treatable with antibiotics. An increasing prevalence of resistance means that to ensure a good treatment outcome it is desirable to test the susceptibility of each infection to different antibiotics. Conventionally, this is done by culturing a clinical sample and then exposing aliquots to a panel of antibiotics, each being present at a pre-determined concentration, thereby determining if the sample isresistant or susceptible to each sample. The minimum inhibitory concentration (MIC) of a drug is the lowestconcentration that inhibits growth and is a more useful quantity but requires each sample to be tested at a range ofconcentrations for each drug. Using 96-well broth micro dilution plates with each well containing a lyophilised pre-determined amount of an antibiotic is a convenient and cost-effective way to measure the MICs of several drugs at once for a clinical sample. Although accurate, this is still an expensive and slow process that requires highly-skilled and experienced laboratory scientists. Here we show that, through the BashTheBug project hosted on the Zooniverse citizen science platform, a crowd of volunteers can reproducibly and accurately determine the MICs for 13 drugs and that simply taking the median or mode of 11-17 independent classifications is sufficient. There is therefore a potential role for crowds to support (but not supplant) the role of experts in antibiotic susceptibility testing.

Tuberculosis is a bacterial respiratory infection that kills about 1.4 million people worldwide each year. While antibiotics can cure the condition, the bacterium responsible for this disease, Mycobacterium tuberculosis, is developing resistance to these treatments. Choosing which antibiotics to use to treat the infection more carefully may help to combat the growing threat of drug-resistant bacteria. One way to find the best choice is to test how an antibiotic affects the growth of M. tuberculosis in the laboratory. To speed up this process, laboratories test multiple drugs simultaneously. They do this by growing bacteria on plates with 96 wells and injecting individual antibiotics in to each well at different concentrations. The Comprehensive Resistance Prediction for Tuberculosis (CRyPTIC) consortium has used this approach to collect and analyse bacteria from over 20,000 tuberculosis patients. An image of the 96-well plate is then captured and the level of bacterial growth in each well is assessed by laboratory scientists. But this work is difficult, time-consuming, and subjective, even for tuberculosis experts. Here, Fowler et al. show that enlisting citizen scientists may help speed up this process and reduce errors that arise from analysing such a large dataset. In April 2017, Fowler et al. launched the project 'BashTheBug' on the Zooniverse citizen science platform where anyone can access and analyse the images from the CRyPTIC consortium. They found that a crowd of inexperienced volunteers were able to consistently and accurately measure the concentration of antibiotics necessary to inhibit the growth of M. tuberculosis. If the concentration is above a pre-defined threshold, the bacteria are considered to be resistant to the treatment. A consensus result could be reached by calculating the median value of the classifications provided by as few as 17 different BashTheBug participants. The work of BashTheBug volunteers has reduced errors in the CRyPTIC project data, which has been used for several other studies. For instance, the World Health Organization (WHO) has also used the data to create a catalogue of genetic mutations associated with antibiotics resistance in M. tuberculosis. Enlisting citizen scientists has accelerated research on tuberculosis and may help with other pressing public health concerns.

The 2021 WHO catalogue of Mycobacterium tuberculosis complex mutations associated with drug resistance: A genotypic analysis.

Background: Molecular diagnostics are considered the most promising route to achieving rapid, universal drug susceptibility testing for Mycobacterium tuberculosiscomplex (MTBC). We aimed to generate a WHO endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction. Methods: A candidate gene approach was used to identify mutations as associated with resistance, or consistent with susceptibility, for 13 WHO endorsed anti-tuberculosis drugs. 38,215 MTBC isolates with paired whole-genome sequencing and phenotypic drug susceptibility testing data were amassed from 45 countries. For each mutation, a contingency table of binary phenotypes and presence or absence of the mutation computed positive predictive value, and Fisher's exact tests generated odds ratios and Benjamini-Hochberg corrected p-values. Mutations were graded as Associated with Resistance if present in at least 5 isolates, if the odds ratio was >1 with a statistically significant corrected p-value, and if the lower bound of the 95% confidence interval on the positive predictive value for phenotypic resistance was >25%. A series of expert rules were applied for final confidence grading of each mutation. Findings: 15,667 associations were computed for 13,211 unique mutations linked to one or more drugs. 1,149/15,667 (7·3%) mutations were classified as associated with phenotypic resistance and 107/15,667 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was >80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were classified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs. Interpretation: This first WHO endorsed catalogue of molecular targets for MTBC drug susceptibility testing provides a global standard for resistance interpretation. Its existence should encourage the implementation of molecular diagnostics by National Tuberculosis Programmes. Funding: UNITAID, Wellcome, MRC, BMGF.

Mycobacterium tuberculosis transmission in Birmingham, UK, 2009-19: An observational study.

Background: Over 10-years of whole-genome sequencing (WGS) of Mycobacterium tuberculosis in Birmingham presents an opportunity to explore epidemiological trends and risk factors for transmission in new detail. Methods: Between 1st January 2009 and 15th June 2019, we obtained the first WGS isolate from every patient resident in a postcode district covered by Birmingham's centralised tuberculosis service. Data on patients' sex, country of birth, social risk-factors, anatomical locus of disease, and strain lineage were collected. Poisson harmonic regression was used to assess seasonal variation in case load and a mixed-effects multivariable Cox proportionate hazards model was used to assess risk factors for a future case arising in clusters defined by a 5 single nucleotide polymorphism (SNP) threshold, and by 12 SNPs in a sensitivity analysis. Findings: 511/1653 (31%) patients were genomically clustered with another. A seasonal variation in diagnoses was observed, peaking in spring, but only among clustered cases. Risk-factors for a future clustered case included UK-birth (aHR=2·03 (95%CI 1·35-3·04), p < 0·001), infectious (pulmonary/laryngeal/miliary) tuberculosis (aHR=3·08 (95%CI 1·98-4·78), p < 0·001), and M. tuberculosis lineage 3 (aHR=1·91 (95%CI 1·03-3·56), p = 0·041) and 4 (aHR=2·27 (95%CI 1·21-4·26), p = 0·011), vs. lineage 1. Similar results pertained to 12 SNP clusters, for which social risk-factors were also significant (aHR 1·72 (95%CI 1·02-2·93), p = 0·044). There was marked heterogeneity in transmission patterns between postcode districts. Interpretation: There is seasonal variation in the diagnosis of genomically clustered, but not non-clustered, cases. Risk factors for clustering include UK-birth, infectious forms of tuberculosis, and infection with lineage 3 or 4. Funding: Wellcome Trust, MRC, UKHSA.

Whole-genome sequencing for surveillance of tuberculosis drug resistance and determination of resistance level in China.

OBJECTIVES: Phenotypic drug susceptibility testing for prediction of tuberculosis (TB) drug resistance is slow and unreliable, limiting individualized therapy and monitoring of national TB data. Our study evaluated whole-genome sequencing (WGS) for its predictive accuracy, use in TB drug-resistance surveillance and ability to quantify the effects of resistance-associated mutations on MICs of anti-TB drugs. METHODS: We used WGS to measure the susceptibility of 4880 isolates to ten anti-TB drugs; for pyrazinamide, we used BACTEC MGIT 960. We determined the accuracy of WGS by comparing the prevalence of drug resistance, measured by WGS, with the true prevalence, determined by phenotypic susceptibility testing. We used the Student-Newman-Keuls test to confirm MIC differences of mutations. RESULTS: Resistance to isoniazid, rifampin and ethambutol was highly accurately predicted with at least 92.92% (95% confidence interval [CI], 88.19-97.65) sensitivity, resistance to pyrazinamide with 50.52% (95% CI, 40.57-60.47) sensitivity, and resistance to six second-line drugs with 85.05% (95% CI, 80.27-89.83) to 96.01% (95% CI, 93.89-98.13) sensitivity. The rpoB S450L, katG S315T and gyrA D94G mutations always confer high-level resistance, while rpoB L430P, rpoB L452P, fabG1 C-15T and embB G406S often confer low-level resistance or sub-epidemiological cutoff (ECOFF) MIC elevation. CONCLUSION: WGS can predict phenotypic susceptibility with high accuracy and could be a valuable tool for drug-resistance surveillance and allow the detection of drug-resistance level; It can be an important approach in TB drug-resistance surveillance and for determining therapeutic schemes.